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ABSTRACT 

The Bullet Cluster is a massive galaxy cluster at z = 0.297 undergoing a major supersonic (Mach 
3) merger event. Using Spitzer 24^m images, IRAC data, optical imaging, and optical spectroscopy, 
we present the global star formation rate (SFR) of this unique cluster. Using a 90% spectroscopically 
complete sample o f 37 non-AGN MIPS confirmed cluster members out to R < 1.7 Mpc, and the 
IRieke et al.l (|2009|) relation to convert from 24/xm flux to SFR, we calculate an integrated obscured 
SFR of 267 M Q yr _1 and a specific star formation rate of 28 M Q yr _1 per 10 14 M Q . The cluster 
mass normalized integrated SFR of the Bullet Cluster is among the highest in a sample of eight other 
clusters and cluster mergers from the literature. Five LIRGs and one ULIRG contribute 30% and 40% 
of the total SFR of the cluster, respectively. To investigate the origin of the elevated specific SFR, 
we compare the infrared luminosity function (IR LF) of the Bullet Cluster to those of Coma (evolved 
to z=0.297) and CL1358+62. The Bullet Cluster IR LF exhibits an excess of sources compared to 
the IR LFs of the other massive clusters. A Schechter function fit of the Bullet Cluster IR LF yields 
L* = 44.68 ± 0.11 ergs s^ 1 , which is ~0.25 and 0.35 dex brighter than L* of evolved Coma and 
CL1358+62, respectively. The elevated IR LF of the Bullet Cluster relative to other clusters can be 
explained if we attribute the "excess" star-forming IR galaxies to a population associated with the 
inf ailing group that have not yet been transformed into quiescent galaxies. In this case, the timescale 
required for quenching star formation in the cluster environment must be longer than the timescale 
since the group's accretion - a few hundred million years. We suggest that "strangulation" is likely 
to be an important process in the evolution of star formation in clusters. 
Subject headings: galaxies: clusters: individual (1E0657-56, The Bullet Cluster) - galaxies: evolution 
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1. INTRODUCTION 

The fraction of star-forming galaxies in galaxy clus- 
ters is known to be suppressed relative to what is typ- 
ically found in the less dense environments of galaxy 
groups and the general field population. The overwhelm- 
ing o bservational evidence (Dressier 1980} iGomez et al.l 
2003) suggests that star formation must be quenched at 
some point during hierarchical mass assembly. Galaxy 
transformations are thought to occur predominantly in 
one of two broad categories - prior to cluster assem- 
bly within the galaxy gro up environment, also known 
as galaxy pre-proce ssing fZabludof f fc Mulchaevi 119981 : 
iKodama et al.l [20011) . or within the cluster environment 
itself. While there is evidence in favor of both scenar- 
ios (e . g. iFuiital 12004} iCortese et aD 120061 : iKovama et"all 
120081: JBerrieret al.1 12009( 1 . in the end the role of the 
cluster environment on the star formation history of its 
galaxies is unclear. Even less clear is the impact of a 
major cluster mer ger on gala xy evolution. 
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fraction of star-forming radio galaxies and a preferential 

1 Department of Astronomy, University of Florida, 
Gainesville, FL 32611-2055; schung@astro.ufl.edu 

2 Department of Physics and Astronomy, Ohio University, 
251B Clippinger Lab, Athens, OH 45701 

3 Alfred P. Sloan Fellow 

4 Harvard-Smithsonian Center for Astrophysics, 60 Garden 
Street, Cambridge, MA 02138 

5 Steward Observatory, University of Arizona, 933 North 
Cherry Avenue, Tucson, AZ 85721 



distribution of emission l ine galaxies between m erging 
subclumps. In contrast, Pogg ianti et al.l (|2004fl found 
that the post-starburst population of dwarf galaxies 
in the Coma cluster lie near the edges of two merging 
substructures, suggesting that the merger quenches 
star formation, though these could be relics of star- 
bursts induced durin g an earlier phase of the merger 
(|Mahaian et al.ll2010l ). 

There has also been much effort placed in under- 
standing the specific physical mechanisms responsible 
for transforming galaxy properties such as morphology 
and star formation. Although there is some observa- 
tional and/or theoretical support for various physical 
processes considered to oc cur in the cluster e nviron- 
ment, su ch as ram pressure dGunn fc Gottlll972|) . stran- 
gulati qn dKawata fc Mulchaevll2008D . and galaxy harass- 
ment (jMoore et al.l 119961 ). it is still unknown which of 
these play a dominant role in transforming a star-forming 
galaxy into a quiescent one. 

This is the second in a series of papers to examine 
the star formation properties of 1E0657-56, also known 
as the Bullet Cluster. The Bullet Cluster is a galaxy 
cluster at z — 0.297 undergoing a major merger event, 
with a collision between the main cluster and subclus- 
ter occuring close to th e plane of the sky with i < 8° 
(|Markevitch et al.ll2004f). A well-defined bow sho ck front 
has been confirmed by iMarkevitch et al.l ([2002D . and is 
propagating through the X-ray gas in the subcluster re- 
gion at a veloci ty of 47 40±630 km s _1 (Mach number 
M = 3.0 ± 0.4) (jMarkevitch 2006). The subcluster itself 
lags behind the shock front due to a wind from the main 
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cluster gas a nd travels at 2700 km s -1 relative to the 
main cluster ([Springel fc Fa rrar 2007]). At this velocity, 
the time since core passage is ^0.25 Gyr ago. 

With our extensive multi-wavelength dataset and op- 
tical spectroscopy, we study star formation out to nearly 
the virial radius of the cluster (R200 = 2.2 Mpc). In 
this series of papers, we will examine various aspects of 
star formation in the Bullet Cluster to understand how 
a recent major merger event affects the star formation 
of cluster galaxies. We also consider the importance of 
different physical mechanisms that may dominate galaxy 
evolution in the diverse environments present in the Bul- 
let Cluster. Star formatio n in the core of the B ullet Clus- 
ter has been studied by iChung et al.l (|2009t ) who have 
found that the strong ram pressure exerted by supersonic 
gas from the cluster merger does not have a significant 
impact on recent star formation in the cluster galaxies. 

In this paper we use the Spitzer Multi band Imag- 
ing P hotometer (MIPS) in the 24/mi band (|Rieke et al.l 
2004), to compare the global star formation rate of the 
Bullet Cluster, as traced by the dusty star-forming galax- 
ies, to those of other merging and non-merging clusters 
from the literature. We also present the infrared lumi- 
nosity function of the Bullet Cluster as a means of in- 
vestigating the nature of the infrared luminous galaxy 
population in this massive major cluster merger. 

2. OBSERVATIONS AND DATA REDUCTION 

The data used in this paper include 3.6/im to 8/im, and 
24/xm imaging from the Spitzer Space Telescope and op- 
tical imaging from the 2.2m ESO telescope. The Spitzer 
data cover a radius of R ~ 1.7 Mpc from the center 
of the Bullet Cluster field (a = 6' l 58 m 26. s 0248, S = 
— 55°56 , 49. 3185), while the optical data extend to R ~ 
4.5 Mpc. We also use optical spectroscopy obtained over 
several years with the Inamori Magellan Areal Camera 
and Spectrograph (IMACS) to determine cluster mem- 
bership for the mid-infrared Spitzer sources based on our 
catalog of 857 known redshifts in the Bullet Cluster field. 
Throughout the paper we assume a cosmology of Hq = 71 
km/s/Mpc, Q M = 0.27, A = 0.73. 

2.1. Spitzer MIPS 

We obtained observations of the Bullet Cluster on 30 
November 2007, in the 24 /im band of the Spitzer Mu lti- 
band Imaging Photometer (MIPS; iRieke et al.|[200l in 
photometry mode. We observed in two cycles, with 14 
frames at 30 seconds per frame for a total integration 
time of ~940 s per pixel in the main exposed region. Us- 
ing a raster map in which frames are dithered along the 
scan and cross-scan directions, with half-array offsets, 
we achieved a total spatial coverage of ~ 12. 6 x 12. 6, 
which corresponds to a radius of 1.7 Mpc from the cen- 
ter of the Bullet Cluster. Individual frames were com- 
bined into a final mosaic with the Spitzer Science Cen- 
ter (SSC) data analys is tool MOsaicker and P oint source 
Extractor (MOPEX; iMakovoz fc Kh an 2005). MOPEX 
also includes a tool for photometry called APEX, which 
performs point-response function (PRF) fitting of point 
sources. Using MOPEX and APEX, we obtained a cat- 
alog of 418 sources detected at a minimum of 5er above 
the background noise. 

We determine the completeness of the MIPS catalog 
via artificial star tests. A bright, unsaturated star in the 



mosaic is used as the point source function (PSF). This 
"PSF" is then duplicated and scaled to create artificial 
stars that cover a range of fluxes. Artificial stars are 
added to the original mosaic, then APEX is used for 
source detection. This process is repeated until a smooth 
completeness function is obtained. The 50% and 80% 
completeness limits are 134/xJy and 310/iJy, with average 
corresponding signal-to-noise (S/N) ratios of ~11 and 
~23, respectively. 

2.2. Spitzer IRAC 

We obtained data from the Sp itzer InfraRed Array 
Camera (IRAC; iFazio et all 12004ft in all four bands - 
3.6/im, 4.5/im, 5. 8/im, and 8.0/im on 14 November 2007 
in full array readout mode with a small scale, cycling 
dither pattern. With eight offset pointings, the total 
coverage is 15' x 15' in the four IRAC bands. 

We created IRAC mosaics from individual frames using 
MOPEX, with source detection and photomet ry done us- 
ing Source Extractor (|Bertin fc Arnoutsl [l996) . We com- 
pute galaxy colors using aperture photometry within a 
6 pixel (3. 6) diameter, which is sufficiently large for ro- 
bust photometry while avoiding contamination from po- 
tential nearby neighbors in a crowded field. No aperture 
corrections were applied because they are minimal when 
considering colors (~0.05 magnitudes). 

2.3. WFI Imaging 

We obtained R, V, and B-band imaging using the Wide 
Field Imager (WFI) on the ESO/MPG 2.2m telescope at 
La Silla, Chile. The data, collected in service-observing 
mode in January 2004, have total exposure times in the 
three filters of 14100s, 6580s, and 5640s, respectively, and 
cover an area of 34' x 34' centered on the Bullet Clus- 
ter field. Individual images have 470s exposure times, 
and were dithered between images to cover the chip 
g aps. We reduced the imag es following the prescription 
of lClowe fc Schneiderl ([2001), with the addition of a pho- 
tometric cor rection based on the scattered light present 
in flat fields ()Koch et al.ll2004T ). Photometric zer o-points 
were determined using standard star fields from Landolt 
(1992). We aligned all images to a common coordinate 
system and construct a photometric catalog using SEx- 
tractor on the R-band image for object detection and 
R-band photometry, and the two-image mode to obtain 
photometry on the B and V-band images. 

2.4. IMACS Spectroscopy 

We have a total of 1122 optical spectra in the Bullet 
Cluster field in the years 2005, 2006, and 2009. In 2005 
and 2006, red sequence galaxies were the primary targets, 
while in 2009 we targetted MIPS sources and blue-cloud 
galaxies beyond the MIPS FOV. We used the Inamori 
Magellan Areal Camera and Spectrograph (IMACS) on 
the 6.5m Magellan Baade telescope with the f/2 cam- 
era. All spectra cover a wavelength range of 4000-9000 
A with a dispersion of 2A per pixel and were reduced 
with COSMOS, the standard data reduction package for 
IMACS spectra. 

To obtain redshifts, we use the IRAF task xcsao to 
cross-correlate spectra with template spectra of four 
galaxy types - giant elliptical, spiral, E+A, and emis- 
sion line galaxy. Of the 1122 spectra, we recover red- 
shifts for 857 galaxies, of which 362 are designated cluster 
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members based on a caustic analysis of the cluster infall 
region (Gonzalez et al. in prep ) , sim ilar to the tech- 
nique employed by Di aferio et al.l ([2005D . To our catalog 
of 362 confirmed m e mbers , we add another 44 members 
from iBarrena et al.1 (j2002h . Figure [1] shows the distribu- 
tion of spectroscopic redshifts, with the cluster members 
highlighted in red. In addition to the Bullet Cluster, 
we find two prominent redshift peaks at z ~ 0.21 and 
z ~ 0.35, containing roughly 90 members each. The Her- 
schel Space Observatory far-infrared and sub-mm prop- 
erties of the background group are studied in relation to 
the Bullet Cluster by Rawle et al (2010, in press). 



3. ANALYSIS 
3.1. MIPS candidate ? 
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To quantify the total star formation rate of the Bullet 
Cluster as traced by 24/xm luminosity, we first identify 
the MIPS sources that are likely to be cluster members. 
We use WFI and IRAC images to exclude galaxies whose 
colors are indicative of either background galaxies or ac- 
tive galactic nuclei (AGN). We then refined our sam- 
ple by spectroscopically targeting MIPS sources that are 
candidate cluster members. 

We start with an initial MIPS catalog of 418 sources 
down to a flux of 43/zJy. We match the WFI (BVR) and 
IRAC catalogs and then cross-match with the MIPS cat- 
alog. Figure [2] shows the R-[4.5] versus B-R colors for 
all MIPS sources with optical and IRAC counterparts. 
Spectroscopically confirmed members from the 2005, 
2006, and first part of 2009 IMACS campaigns demon- 
strated that the star-forming cluster members form a 
tight diagonal locus in the R-[4.5] vs B-R space. We de- 
fine our cluster candidate sample as galaxies within the 
two solid diagonal lines shown in Figure [2j and blueward 
of B-R=3. 

We then use IRAC colors to i dentify AGN via the 
"AGN wedge" (ILacv et alJl200i IStern et alJl2005h . as 
shown in Figure [3j We exclude all galaxies within the 
AGN wedge from further analyses because their mid- 
infrared luminosity may be dominated by AGN activity 
rather than by star formation. We also exclude X-ray 
AGN using a catalog of 145 X-ray point sources extracted 
from Chandra data that cover the central ~ 20' x 20', 
overlapping the entire MIPS FOV. The catalog contains 
X-ray sources down to a flux of 2.5 x 10~ 16 ergs cm~ 2 
s^ 1 in the 0.5-2 keV band, which is a luminosity of 
L x = 7x 10 40 ergs s _1 at the redshift of the Bullet Clus- 
ter. Among the MIPS confirmed cluster members and 
candidates, we exclude all X-ray point sources with an 
X-ray luminosity of L x > 10 41 ergs s _1 . 

In addition to IRAC and X-ray selection of 
AGN, we also utilize our optical spectroscopy to 
construct Baldwin-P hillips- Terlevich diagrams (BPT; 
iBaldwin et all 119810 as illustrated in Figure H The 
dotted and dashed lines indicate the boundaries from 
iKewlev et al.l (|2006[ ) that separate purely star forming 
galaxies from Seyferts and LINERs. Note that many of 
the MIPS sources are missing in Figure |4] because one 
or more of the four required emission lines could not be 
measured, often due to a prominent sky emission line 
that appears at the same wavelength as H/3 at the Bullet 
Cluster redshift. 

Figure H] reveals only two MIPS sources classified as a 



Seyfert or LINER that have not been identified as AGN 
using either IRAC colors or X-ray emission. One of these 
is a ULIRG (star symbol), which appears close to the 
Seyfert/LINER boundary when using the [SII]/Ha ratio, 
and is classified as a LINER using the [OI] /Ha ratio. 

The two confirmed MIPS sources in the IRAC AGN 
wedge (Figure [3]) are not shown on the BPT diagram 
because they lack the necessary emission lines. However, 
of the 10 confirmed IRAC AGN wedge sources that are 
not MIPS members (or beyond the MIPS FOV), two are 
plotted in Figure |4j with one classified as a LINER, and 
the other as an HII dominated galaxy. 

Of the three X-ray point sources in Figure |31 one is 
illustrated in Figure|l]as a large square and appears close 
to the boundary of HII dominated galaxies and Seyferts. 

In total, there are 8 AGN identified among the con- 
firmed MIPS sample, using the three methods - IRAC 
colors, X-ray emission, and optical emission line ratios. 
Of these eight, two are identified solely from the BPT 
diagnostic, one classified as a LINER and the other a 
Seyfert. The LINER is included in our sample, assuming 
that 65% of its IR flux is powered by star formation (see 
§4.ip . The Seyfert, which is excluded from our sample, 
would contribute a neglible fraction to the global SFR. 

In addition to isolating the AGN population, Figure [3] 
also illustrates that our selection of cluster candidates 
based on R-[4.5] and B-R colors is an effective way to 
cull interlopers. The left and right panels of Figure [3] 
show the IRAC color distribution before and after we ap- 
ply the R-[4.5] and B-R color selection, respectively. We 
show the model color s of M82 (a local starburst galaxy) 
at z= 0,0.3,0.5, and 1 (|Devriendt et al.lll999HStern et al.l 
2005) to illustrate that our color selection is effective in 
removing galaxies whose IRAC colors are consistent with 
those of a high redshift starburst (open red circles near 
z ~ 1 to 2). 

The right panel of Figure [3] highlights the AGN sources 
among the MIPS sample, including four galaxies in the 
IRAC AGN wedge, three X-ray point sources outside the 
AGN wedge, and one Seyfert identified from the BPT di- 
agram, which is semi-hidden behind the ULIRG symbol. 
Among the non-AGN MIPS cluster members, there are 
two outliers, both of which do not fit in the diagonal 
star forming sequence (roughly outlined by the LIRGs 
and ULIRG shown as star symbols) nor the locus of 
passive early type galaxies near [3.6]-[4.5]~0.2 and [5.8]- 
[8.0]~0.2. One outlier has a color of [3.6]-[4.5]< 0, due to 
a blending of sources in the 3.6/im and 4.5/im bands. It 
is unclear why the second outlier has IRAC colors consis- 
tent with those of a high redshift starburst. Perhaps an 
AGN component undetected in X-ray is contaminating 
its IRAC colors. We include both sources in the following 
analyses since they are cluster members with 24/xm emis- 
sion and have no sign of a dominant AGN component. 
Both galaxies have a SFR< 3.6 M©yr _1 and contribute 
negligibly to any results presented in this paper. 

3.2. Spectroscopic Completeness 

We match our photometric catalogs to our spectro- 
scopic one, which includes 857 redshifts within the full 
IMACS FOV (25' x 25'), centered on the Bullet Cluster 
(Figure Q}. Of these, 362 galaxies are confirmed clus- 
ter members and 495 are interlopers. Combined with 
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the iBarrena et al.l (|2002D catalog, we have a total of 406 
confirmed cluster members. 

In estimating the spectroscopic completeness of MIPS 
galaxies, we consider only the 139 MIPS sources that lie 
within the diagonal color boundaries in Figure [2] Out of 
those, nine lie within the AGN wedge - two confirmed 
members, five interlopers, and two cluster candidates. 
The spectroscopic completeness within the AGN wedge 
for MIPS sources is -80%. 

Out of the 130 MIPS galaxies that lie outside the AGN 
wedge, we have redshifts for 118 galaxies, with 39 con- 
firmed cluster members, 78 interlopers, and 12 cluster 
candidates with no spectroscopy. We are 90% spectro- 
scopically complete for MIPS sources outside the AGN 
wedge, with —35% being cluster members. 

There are a total of 43 confirmed members with MIPS 
emission, including two cluster members whose R-[4.5] 
colors are slightly beyond the color selecti on of Figure [21 
and fiv e members whose redshifts are from lBarrena et al.l 
(2002). Among these 43, we exclude two in the IRAC 
AGN wedge, three that are X-ray point sources, and one 
that is a Seyfert. Of the additional 14 MIPS cluster can- 
didates chosen based on IRAC and optical colors, two 
are in the AGN wedge and hence excluded. Our final 
non-AGN MIPS sample consists of 37 confirmed mem- 
bers and 12 candidate members, of which we expect —4 
to be members based on our success rate for targeting 
MIPS cluster members. 

3.3. Total Infrared Luminosity 

To determine the total infrared luminosity Ljr (A = 
8 — 1000 /jm) and star formation rate (SFR) from the 
observed 24 / im flu x /(24), we use the prescription in 
iRieke et al.l ((2009), who construct model spectral en- 
ergy diagrams (SEDs) for a sample of local LIRGs and 
ULIRGs and derive relations between 24/im flux, star for- 
mation rate, and total infrared luminosity. The scatter 
a ssociated with the L 24 — Ljr relation is 0.13 dex. 

iRawle et al.l ((2010h have found that for a sample of 
23 Bullet Cluster galaxies, —30% ar e found to have an 
exc ess of 100/xm flux rel ative to the IRieke et al.1 (2009) 
and lDale fc Heloul (|2002t ) templates. Although the devi- 
ation between Herschel derived and 24/im derived total 
infrared luminosities can be up to a factor of 4 for some 
sources, these are typically low-luminosity sources fainter 
than our S/N criterion for the 24 micron sample. Among 
the galaxies used in this paper (with available Herschel 
data), the ratio of SFR^/fl to SFR24 is close to 1, with 
a scatter of 0.6, indicating that the 24/im flux yields ac- 
curate total IR luminosities and star formation rates for 
our sample. 

Although we have a decent understanding of the 
systematic u ncerta inty in the Lm derived from the 
IRieke et al.1 ((2009) templates, much of the previous 
work in the l iterat ure uses the galaxy templates of 
iDale fc Heloul (|2002l) to derive Lm from observed 24/im 
flux. In the left panel of Figure we show the rela- 
tion^ betweeD_JLj_B and ob served 24/im flux derived from 
the IDale fc Heloul (2002) templates at z=0.3 for our 49 
non-AGN MIPS members and c luster candidates. We 
compare the lDale k Heloul (120021) Lib, values to those de- 
rived using the IRieke et al.l ([2009) calibration. The lower 
portion of the left panel in Figure [5] shows the difference 
between the two sets of Lm as a function of observed 



24/im flux. 

A mong the 49 MIPS sources, the IRieke et all (|2009fi 
and Dale fc Heloul (|2002l ) templates yield a difference in 
total infrared luminosity of approximately -0.5 dex, -0.25 
dex, and 0.1 dex for the low, median, and high end of 
24/im fluxes in our sample, respectively. The brightest 
24/im source shown in Figure [5] is 65% of the ULIRG 
flux. 

3.4. Star Formation Rate 

We calculate the star formation ra te from the observe d 
24/im flux by using equation 14 in IRieke et al.l (2009) . 
The uncertainty in the SFR derived from this relation 
is —0.2 dex. This error is dominated by the scatter in 
the L24 — Ljb. relation, with an additional error due to 
the underlying assumption that most of the young stel- 
lar light is absorbed and re-radiated in the infrared. Al- 
though an average correction has been applied to account 
for the loss of young stellar light directly to the ultravi- 
olet (and therefore never seen in the infrared) , there is a 
variation from galaxy to galaxy that adds an additional 
uncertainty to the star formation rate calibration. The 
0.2 dex estimate does not include uncertainties in the 
assumed initial mass function or the theoretical relation 
between SFR and Lir. We assign an error of 0.2 dex to 
all star formation rates calculated in this paper. 

In addition to the uncertainty in the star formation 
rate, we need to understand any sy stematic offs e ts be - 
tween the SFRs calculated from the IRieke et"aT1 (j2009[ ) 
relations, versus the SFRs quoted in much of the pre- 
vious literature using the IDale fc Heloul (120021) galaxy 
templates and the original TKennicuttl (|1998|) SFR-Ltr 
relation. Star formation r ates obtained from the SFR- 
Ljr relation presented in IRieke etldl (|2009l ) are sys- 
tematically lowe r by —0.2 dex relative to the standard 
lKennicutn(|1998| ) relation due to a difference in the initial 
mass function (IMF). An unbroken Salpeter I MF of slope 
-1.35 f rom 0.1 to 100 M Q is used in the original|Kenmcutt] 
(|1998f ) derivation, whereas more recent studies suggest 
that a Salpeter-like IMF with a shallower slope for the 
low mass end is more applicable for extragalactic star- 
formi ng regions (jRieke et al.lll993t lAlonso-Herrero et al.1 
120011 ). This broken Salpeter IMF has a total mass that is 
—0.66 times the mass of a single power-law Salpe ter IMF, 
and y ields a similar proportion of s tars to the iKroupal 
(pOOa) IMF and the lChabrierl (l2003h IMF. 

The right panel of Fig ure [5] shows the co mparison of 
SFRs derived from the IRieke et al.l (12009)) calib ration 
compared to those from th e IDale fc Heloul (|2002D tem- 
plates and lKennicutti (|1998t ) relation, as a function of the 
observed 24/im flux of the Bullet Cluster MIPS sources. 
The SFRs calculated from the two methods are systemat- 
ically offset from each other more significantly than in the 
case of total infrared luminosities. The bot tom portion 
of the right panel of Figure [5] shows that the IRieke etldl 
(2009) calibration yields systematically lower star for- 
mation rates by a factor of —2.6 to 1.2 for the range of 
24/im fluxes in the Bullet Cluster sample, as noted by 
the dotted vertical lines. The difference in SFRs derived 
from the two methods is in part due to the systemat- 
ically lower infrar ed luminosities derived from the the 
IRieke etldl (|2009| ) templates, as well as t he as sumption 
of a different IMF from the iKennicuttl (J1998D relation 
which produces fewer low mass stars. 
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4. RESULTS & DISCUSSION 
4.1. LIRGs and ULIRG 

Among the 37 non-AGN MIPS confirmed cluster mem- 
bers, five are LIRGs (L/_r > 10 11 L Q ) and one is a 
ULIRG (L IR > 10 12 L Q ). Although the ULIRG docs 
not appear within the AGN wedge nor as an X-ray point 
source, its optical spectrum exhibits broad emission lines 
and emission line r atios on the BPT diagrams (BPT; 
I Baldwin et al .1119811 ) that classify this galaxy as a border- 
line Seyfert/LINER when using the [SIIJ/Hq ratio and a 
LINER with the [01] /Ha ratio (Figure g]). 

Mid-infrared and far-infrared emission in ULIRGs 
arises from dust that is heated by a young stellar pop- 
ulation and/or a central AGN. In the case of LINERs, 
some studies have found that the infrared emission is 
dominate d by star f ormatio n rather than nuclea r activ- 
ity (e.g. [Saty^SiEinilOOll^nii^iEinilOOi). Using 
six different method s to determine AGN contribution, 
iVeilleux et~aT1 (|2009f) found that the AGN contribution 
to the bolom etric luminosi t y in L INERs is ^15% to 35%. 
In addition, ISturm et all (J2006D found that IR-bright 
LINERs tend to have mid-IR SEDs consistent with star- 
bursts, whereas SEDs of IR-faint LINERs appear more 
AGN dominated. The LINER in our sample has a total 
IR luminosity of ~ 10 12 L , making it more likely to be 
starburst dominated. While we do not know the precise 
AGN contribution in our ULIRG, we adopt a value of 
35%, in a cautious effort not to overestimate the true 
star formation rate. In the proceeding analyses and dis- 
cussion, we include this ULIRG in our sample, but always 
with a star formation rate that assumes that only 65% of 
the infrared flux is powered by star formation. The sky 
coordinates of the five LIRGS and one ULIRG are listed 
in Table [H along with their total infrared luminosities 
and star formation rates. 

Figure |5] shows the specific star formation rate (SSFR) 
and stellar mass for the LIRG and ULIRG population, 
and the remaining non-AGN MIPS sample of 37 con- 
firmed members and 12 cluster candidates. There is 
a clear trend of decreasing specific SFR with increas- 
ing stellar mass. Although there is an inherent selec- 
tion bias against low mass, low specific SFR galaxies, 
there is evidence that the general correlation seen in Fig- 
ure is real. The upper envelope of galaxies with high 
SFRs (parallel to and above the 50% completeness line in 
Figure [5J) should represent a nearly complete sample of 
bright 24/^m sources. Among this nearly complete sam- 
ple of high SFR galaxies, the trend of decreasing specific 
SFR as a function of stellar mass is still apparent. In 
addition, the correlation seen in Figure [6] has been ob- 
served in several other studies, in both field and cluster 
galaxi es (Fcul ner et al.ll2005HPerez-Gonzalez et alj|2005t 
INoeske et al7 T2007; Vul cani et all 120101) . We note that 
Figure H] is intended to illustrate a general trend between 
SSFR and stellar mass (not a precise fit or relation), as 
well as the range of SSFR and stellar mass values in our 
MIPS sample. 

The five LIRGs and ULIRG, which are highlighted 
(star symbols) in Figure [51 represent an outlier popula- 
tion among the MIPS sources. They show excess levels of 
star formation given their stellar mass, with the ULIRG 
being the most extreme outlier from the nominal SSFR 
versus stellar mass relation. 



All galaxies listed in Table[2]lie beyond R~l Mpc from 
the cluster center, with the two most luminous galaxies 
located at R~1.7 Mpc. The intermediate to cluster out- 
skirts region is an environment similar to that of galaxy 
groups, where galaxy-galaxy interactions are more likely 
to occur than in high density cluster cores or in the low 
density field. However, the visual morphologies of the 
LIRGs do not show any obvious signs of recent merger 
activity, based on BVR imaging from the Magellan 6.5m 
Baade telescope. At least two out of the five LIRGs are 
clearly spiral galaxies, while morphologies of the other 
three are not entirely clear, though they do show some 
faint spiral structure. Figure [7] shows an image of the 
ULIRG taken with the Advanced Camera for Surveys 
on the Hubble Space Telescope in the F606W filter (P.I. 
Holland Ford; Proposal ID 10996). It is a barred spi- 
ral galaxy, with no obvious signs of recent interaction or 
mergers, though there is some assymetry in the spiral 
arms. 

All six LIRGs/ULIRG have optical colors consistent 
with those of star-forming, blue cloud galaxies, as shown 
in the color magnitude diagram in Figure [5] The IRAC 
color-color diagram (Figure[3]) also shows that the LIRGS 
and ULIRG have colors consistent with star-forming 
galaxies, with the ULIRG having a similar color as the 
starburst galaxy M82 evolved to z = 0.3. 

Table [2] also lists the six cluster members that are clas- 
sified as AGN from IRAC, X-ray, and optical spectro- 
scopic data. One out of the six galaxies is an X-ray point 
source and LIRG, with L[r ~ 10 11 L©, lying close to the 
IRAC AGN wedge. Optical colors of the six confirmed 
and two candidate member AGN are shown in Figure |SJ 

4.2. Global Star Formation Rate 

The integrated SFR of the 37 MIPS confirmed cluster 
members within R < 1.7 Mpc is 267 M Q yr~ 1 , with a 0.2 
dex uncertainty (see H3.4I) . The 12 MIPS candidate mem- 
bers contribute an additional 27 M Q yr _1 . Normalizing 
by the cluster mass M(< 1.7 Mpc) = 9.5 x 10 14 M Q , de- 
rived from the caustic analysis of the cluster infall region 
(Gonzalez et al in prep), we obtain 28 M Q yr _1 per 10 14 
M Q for the confirmed members sample, with the ULIRG 
contributing 40% and the five LIRGs an additional 30% 
(see Table©. 

The Bullet Cluster integrated star formation rate nor- 
malized by cluster mass is presented in Figure [HI 
along with known values for eight other clusters taken 
from IGeach et all (|2006[) and references therein, and 
IHaines et all (|2009D . who use 24^m MIPS or 15/im ISO- 
CAM data to calculate obscured SFRs. In order to fa- 
cilitate a fair comparison of SFR in the Bullet Cluster 
and the other clusters shown in Figure |H1 we impose an 
infrared luminosity lower limit on our MIPS sample. In 
addition, we scale the st ar formation rate t o account for a 
system atic offset in the iRieke et all (|2009l ) relations ver- 
sus the IDale fe Hel ou (2002) templates with the original 
IKennicutti (119981) rela tion (Figure EJ. 

IGeach et all (2006) apply a lower limit of Lju = 
6 x 10 10 L Q (SFR ~ 10 Moyr" 1 ), for all their clus- 
ters with the exception of MS0451-03. This lower limit 
was calculated for a 2 00 /zJy source at z = 0. 39 using the 
IDale fc Heloul (|2002f) templates. Using the IRieke et all 
(2009) calibration, the same source at z = 0.39 yields 
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a total IR luminosity of L IR = 3.8 x 10 10 L (SFR = 
3.6 MQyr -1 ). After applying a lower limit of Lj R = 
3.8 x 10 10 IjQto our full sample of 37 confirmed members 
and 12 cluster candidates, 20 galaxies remain - 19 mem- 
bers and one candidate. We then scale the SFR of each 
galaxy using the relation presented in the bottom right 
panel of Figure \Sj in orde r to convert S FRs calculated 
from t he lRieke et al.l (12009ft relation t o thelDale fc Heloul 
([2002ft system used by IGeach et~aT1 ([20061 ). After mak- 
ing these adjustments, we obtain an integrated SFR of 
445 M Q yr~ 1 for the 20 galaxies, and a cluster mass nor- 
malized SFR of 47 M Q yr _1 per 10 14 M Q , as shown in 
Figure [SJ 

All star formation rates illustrated in Figure [9] are de- 
rived from mid-infrared data, primarily with the MIPS 
24fim band and in some cases supplemented with data 
from the Infrared Satellite Obser vatory (ISO) 15jjm 
band. Integrated SFRs calculated bv IGeach et al.l (J2006I ) 
include mid-IR sources out to a cluster radius of R < 2 
Mpc, while the SFR of A1758 is obtain ed fro m mid-IR 
data out to R < 3 Mpc ([Haines et al.l ((2009)). This is 
comparable to the R < 1.7 Mpc area that we survey. 

The mass normalized SFR of the Bullet Cluster pre- 
sented in this paper is a more precise estimate of the SFR 
than what was possible for many of the clusters shown in 
Figure [9j which mostly lacked mid- infrared data beyond 
the cluster core. Extrapolation of mid-infrared sources 
and field subtraction were applied in such cases, to place 
upper and lower limits on the global specific SFR, illus- 
trated as horizontal bars and arrows on Figure [9] In 
addition, the spectroscopy was sparse for several of the 
clusters presented in Figure H3 whereas the Bullet Cluster 
SFR is based on a 90% spectroscopically complete MIPS 
sample. 

While the integrated SFR of the Bullet Cluster ex- 
cludes contamination from both mid-IR and X-ray 
AGN man y of the clusters co nsidered by IGeach et al.l 
(2006) and lHaines et al.l (|2009[) lack the necessary data 
for a simila r ly co mplete AGN removal. In A1758, 
IHaines et~ al. (2009) identify two X-ray point sources that 
are removed as AGN, though no consideratio n is made for 
mid-in frared AGN without X-ray emission. IGeach et al.l 
(2006J) remove AGN from CL0024+16 and MS0451-03 by 
excluding mid-IR sources that have optical and K-band 
colors similar to the e xpected colors of pas sive early type 
galaxies (E/S0) from King fc Ellid (fl985ft models. Al- 
though this may help eliminate AGN with early type 
hosts, there are two main issues with using this color 
criterion to select out AGN. 

The first problem is that identifying AGN solely with 
this method can be significantly incomplete. As our own 
data show in Figure [31 none of the AGN (six confirmed 
members, two cluster candidates) have IRAC colors con- 
sistent with the locus of early type galaxies near [3.6]- 
[4.5]~0.2 and [5.8]-[8.0]~0.2." In Figure M the optical 
color-magnitude diagram shows that only two out of the 
eight AGN are on the red sequence. 

The second problem is that attempting to select out 
early-type AGN hosts based on their optical and K-band 
colors can inadvertently remove many dusty star-forming 
galaxies, which may actually be spiral galaxies with op- 
tically red colors. Figure [5] shows that eight out of the 37 
non-AGN MIPS cluster members lie on the red sequence. 



Although 20% of the non-AGN MIPS cluster members 
are optically red, these galaxies have a total SFR of ^14 
M Q yr _1 , contributing only ~5% to the integrated SFR. 
While the global SFR in the Bullet Cluster is dominated 
by LIRGs/ULIRG, optically red mid-IR galaxies may 
have more of an impact on the global SFR in clusters 
that lack a LIRG/ULIRG population. In general, it is 
important not to exclude galaxies as AGN based purely 
on their optical colors, as we have seen that up to 20% 

of star-forming galaxies can lie on t he red sequen ce. 

Had we used a similar criterion as IGeach et al.l {2006) 
to exclude AGN based on IRAC and optical colors of pas- 
sive galaxies, we would not have detected any of the six 
confirmed or two candidate AGN, whose total infrared lu- 
minosity converts to a star formation rate of 42 M Q yr _1 
for the confirmed members, and 8 M^yr" 1 for the two 
cluster candidates, using the iRieke et al.l (2009) calibra- 
tion. Adjust i ng for the different SFR calibration used in 
IGeach et al] (|2006l) . the global specific SFR of the Bullet 
Cluster would increase from 47 to 57 M Q yr _1 per 10 14 
M Q in Figure [9] 

4.3. Specific SFR in the Bullet Cluster 

The mass normalized integrated SFR of the Bullet 
Cluster is among the highest in the sample of local to 
intermediate redshift clusters shown in Figure [9j second 
only to CL0024+16. The specific SFR of the Bullet Clus- 
ter is also comparable to that of A1758, a cluster merger 
at nearly the same redshift. Interestingly, the three most 
active clusters (CL0024+16, A1758, and the Bullet Clus- 
ter) in Figure|9]are all cluster mergers. However, all three 
arc in different merger states, with different masses and 
dynamical histories. Although Figure|H]is somewhat sug- 
gestive that cluster mergers have elevated global specific 
SFRs compared to non-merging clusters, a larger sample 
of clusters is necessary to confi rm a trend. 

Previous studies from iKodama et al.1 (12 004); 
IGeach et all (|2006f l: iBai et al.l ([2007ft : iKovama et al.1 
(2010) have suggested that there is a redshift dependence 
on the global star formation rate of clusters, anywhere 
from (1 + z) A to (1 + z) 7 , based on plots similar to 
Figure [HI Although Figure [H shows some suggestion of 
evolution in the specific SFRs, there is clearly much 
scatter between individual clusters. We emphasize that 
the integrated SFR is highly sensitive to small number 
statistics. In the case of the Bullet Cluster, one galaxy 
is contributing ^40% of the entire integrated SFR. 

Another way to compare the star formation activity in 
different clusters is to look at the mass normalized inte- 
grated SFR as a function of cluster m ass. As found in 
IBai et al.l (|2007ft ; IKovama et al.l ([2010ft . there is a depen- 
dence of specific SFR on total cluster mass, although it 
does not display a much stronger correlation than with 
redshift. 

Of all the clusters in Figure H CL0024+16 and the 
Bullet Cluster have the highest number of spectroscop- 
ically confirmed 24/mi MIPS sources, with 45 and 44 
cluster members respectively, prior to excluding AGN. 
Although CL0024+16 has a specific SFR that is a fac- 
tor of ~4 higher than the Bullet Cluster, this could 
in part be driven by its small cluster mass (M(< 
2Mpc) = 6.1 x 10 14 ; iKneib et al.1 12003ft . as several au- 
thors have noted a correlation between global specific 
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SFR and cluster mass (e .g. iFinn et all 120051 iBai et al.1 
[20071 iKovama et al .112010ft . 

Among the three clusters with highest specific SFR 
(A1758, Bullet Cluster, and CL0024+16), we note that 
the correction for AGN contamination in CL0024+16 
and A1758 is less complete than for the Bullet Cluster, as 
is the spectroscopic completeness of A1758. In addition, 
the SFR for A1758 is integrated from candidate members 
out to R < 3 Mpc, compared to R < 1.7 Mpc for the 
Bullet Cluster. It has been shown in several studies that 
infrared and submm luminous galaxies with high levels of 
star formation rate are preferential l y distributed in clus - 
ter outskirts (e.g. iMa et all 12015 feraglia et al.l 120101 ) . 
Therefore a comparison of integrated SFR (not normal- 
ized by mass), within R < 3 Mpc versus R < 1.7 Mpc is 
not necessarily a fair view on the relative star formation 
activity between two clusters. 

The integrated SFR of the Bullet Cluster is strongly 
driven by a small number of infrared luminous galaxies. 
Among our sample of 37 MIPS confirmed members ex- 
cluding AGN, one ULIRG contributes 40% of the total 
SFR, while five LIRGs contribute an additional 30%. We 
will see in the subsequent sections that these few galaxies 
not only dominate the global SFR, but have a significant 
impact on the infrared luminosity function as well. 

4.4. Infrared Luminosity Function 

To further investigate the origin of the Bullet Cluster 
elevated global star formation rate relative to other clus- 
ters, we construct an infrared luminosity function of the 
star-forming population and compare it to the IR LF of 
non-merging clusters with similar redshift and mass. 

The infrared luminosity function of the Bullet Cluster 
is presented in Figure 1 101 with la Poisson error bars. 
We include two versions of the Bullet Cluster IR LF - 
one which excludes all known AGN wedge and X-ray 
point sources (left panel), and one which does not (right 
panel) . The aim of the latter IR LF is to illustrate the 
impact of AGN contamination. The first sample contains 
39 confirmed cluster members and 12 cluster candidates, 
while the latter sample has 44 confirmed members and 14 
cluster candidates. In both cases, we weight the cluster 
candidates by 0.35, which is the probability that a MIPS 
source within the color boundaries set in Figure [5] is an 
actual cluster member. We apply a completeness cor- 
rection derived from artificial star tests for data fainter 
than the MIPS 80% completeness limit, shown as a dot- 
ted vertical line in Figure fTUl 

Overplotted on FigureHUJare the IR LFs of CL1358+62 
(z = 0.328) a n d the Co ma cluster (z = 0.024) 
from iTran et all (|2009l ) and IBai et al.l (|2009D . Both 
CL1358+62 and Coma have approximately the same to- 
tal mass as the Bullet Clust er (M~ 10 15 M g : lKubo et al.1 
[20071 Hoekstr aet al.lll998ft . The Coma IR LF shown 
in Figure QJJ] is evolved to z = 0.3 using the relations 
L* H oc (1 + z) 32 and <S>* IR ex (1 + z) 17 from IBai et "all 
(2009). In addition to the CL1358+62 and Co ma, we also 
overpl ot the IR LF of SG1120 {z = 0.37) from lTran et~aTI 
(2009). SG1120 consists of four galaxy groups in the pro- 
cess of cluster as sembly, predicted to form a Coma-like 
cluster by z = ([Gonzalez et al.ll2005l) . 

The infrared luminosity function of Coma and 
CL1358+62 are derived from 24^m MIPS data that cover 
^3x3 and ~ 2.5 x 2.5 Mpc, respectively. This is compa- 



rable to the Bullet Cluster MIPS data, which covers an 
area of - 3.4 x 3.4 Mpc. The spatial coverage of SG1120 
is ~ 6 x 6 Mpc because the four galaxy groups require a 
more extensive survey area than the clusters. In addition 
to the similar spatial coverage of CL1358+62, Coma, and 
the Bullet Cluster, we emphasize that the total mass of 
all systems shown in Figure ITU] are comparabl e to the 
mass of the Bullet Cluster. ITran et al.1 (|2009l ) do not 
take AGN contamination into a ccount in the IR LFs of 
CL1358+62 and SG1120, while IBai et~aTl ([20091 ) cross- 
matched their Coma data with the Catalogue of Quasar s 
and Active Galactic Nuclei ([Veron-Cettv fc Veronl 2003). 

We fit a Schechter function to the Bullet Cluster data, 
fixing the faint en d slope to a = ~ 1-4, adopted from the 
IR LF of Coma ([Bai et al. 1 12009ft . The ULIRG, which 
occupies the brightest bin in Figure [T0l is excluded from 
the fit because it is a clear outlier from the smooth dis- 
tribution of IR star-forming galaxies. 

Results from fitting a Schechter function to the Bullet 
Cluster data with and without AGN are summarized in 
Table [3[ along with the Schechter parameters for Coma, 
CL1358+62 and SG1120. We note that the Schechter pa- 
rameters for the Bullet Cluster remain consistent within 
lcr whether or not we allow a to vary during the fit. 

Figure [TU] shows that the Bullet Cluster luminosity 
function is enhanced relative to the IR LFs of CL1358+62 
and Coma at all values of Lir. In addition to showing 
a higher surface density of IR star-forming galaxies in 
all Lin bins, the IR LF of the Bullet Cluster extends to 
log Lib, ^ 45 (not including the ULIRG), whereas the IR 
LFs of CL1358+62 and Coma truncate at \ogL IR ~ 44.4 
and log Lir ~ 44.6, respectively. For the IR LF ex- 
cluding all known AGN, we obtain L* = 44.68 ± 0.11,, 
which is ^0.25 and 0.35 dex brighter than the L* value 
of evolved Coma and CL 1358+62, respectively. 

Figure Q7J] also illustrates that inclusion of even a small 
number of AGN sources can have a noticeable impact on 
the IR LF, particularly at the bright end. If we include 
all known AGN, (six confirmed members and two clus- 
ter candidates weighted by 0.35 - the probability that 
they are cluster members), L* increases to 44.81 ± 0.13, 
and is ^0.4 and 0.5 dex brighter than L* of Coma and 
CL1358+62, respectively. 

Although the Bullet Cluster IR LF has a value of 
L* that is brighter than that of evolved Coma or 
CL1358+62, it is still fainter than L* of SG1120 by 
0.30 (0.18) dex, excluding (including) the AGN popula- 
tio n in the Bu l let Cl uster sample. It has been shown 
by ITran et al.l (|2009f) that SG1120 exhibits an excess 
of 24/zm sources compared to CL1358+62 and evolved 
z ~ clusters, indicating that the infrared galaxies of 
SGI 120 represent a galaxy group population whose star 
formation has not been quenched to typical cluster lev- 
els. We see from Figure ITOl that while the Bullet Cluster 
has a higher number density of IR star-forming galaxies 
compared to other clusters, it is still suppressed relative 
to these galaxy group IR LFs. 

One explanation for the excess of IR bright sources in 
the Bullet Cluster relative to other clusters can be the 
presence of galaxies belonging to the infalling group pop- 
ulation associated with the bullet, which may contribute 
significantly to the overall elevated IR LF. In this case, 
the IR LF would be representative of a combined galaxy 
group and galaxy cluster environment. To test this see- 
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nario, we scale down the SG1120 IR LF by the subcluster 
to SG1120 mass ratio, thereby creating a rough proxy for 
a galaxy group IR LF at subcluster mass. 

To determine the subcluster mass, we assume 
a main cluster to su bcluster mass ratio of 1:10 
(jSpringel fc Farrarl I2007J ) and a total cluster mass of 
M 2 oo = 1.5 x 10 15 (Gonzalez et al. in prep). This 
yields a subcluster mass of M sub = 1.4 x 10 14 M©. With 
a lower limit on the mass of SG1120 (5.3 x 10 14 M Q ; 
IGonzalez et alJl2005l ). the subcluster to SG1120 mass ra- 
tio is approximately 1:4. Figure [Til shows the IR LF of 
the Bullet Cluster overplotted with the IR LF of evolved 
Coma added with the IR LF of SG1120 after it has been 
scaled down by a factor of four. Overall, the IR LF of 
combined Coma and scaled SG1120 is a good match to 
the observed Bullet Cluster IR LF, with the two LFs hav- 
ing a similar effective L* . This supports the idea that the 
Bullet Cluster IR LF is a combined distribution, repre- 
senting galaxies from both a cluster and group popula- 
tion. 

Quantitatively, the Bullet Cluster LF is still a factor 
of ~1.5 higher than the evolved Coma+scaled SG1120 
model at most luminosity bins. There are several poss- 
sible explanations. First, the intrinsic variation in the 
IR LF between groups remains poorly constrained-one 
would only need a factor of ~1.5 higher surface density 
in the subcluster to reproduce the observed Bullet Clus- 
ter IR LF. Second, the excess may be evidence of trig- 
gered star formation in the Bullet Cluster due to the 
merger. In this case the factor of 1.5 can be considered 
an upper bo und on the t o tal en hancement induced by 
the merger. iChung et all (|2009f ) found that ram pres- 
sure associated with the supersonic shock front does not 
have a significant impact on recent specific star forma- 
tion rates of galaxies within the central R^ 0.5 Mpc re- 
gion. However, this does not exclude the possibility that 
galaxies in the cluster outskirts may experience low level 
triggering of star formation due to ram pressure effects 
upon first entering the cluster merger environment. Since 
ram pressure is a fast-acting mechanism that can lead to 
quenched star formation in as short a time as ~100 Myr, 
galaxies with triggered star formation in the outskirts re- 
gion could already be transformed into quiescent galaxies 
with little to no star formation, by the time they reach 
the cluster core. 

In cither case, the overall enhancement of star- 
formation in the Bullet Cluster relative to quiescent clus- 
ters of similar mass appears to be dominated by the pop- 
ulation associated with the infalling galaxy group. This 
is supported by the two IR LFs shown in Figure[TTJ which 
are consistent within a factor of 1.5. Group to group 
variation in IR LFs and/or a small effect from the clus- 
ter merger itself may also contribute a small excess of IR 
star-forming galaxies in the Bullet Cluster. 

4.5. Lessons Learned 

We have seen throughout the previous sections that 
although the integrated specific SFR can be a good mea- 
sure of the overall level of activity in a cluster, we must 
exercise caution when comparing this quantity among 
different clusters. As described in N4.21 the contamina- 
tion from an unknown AGN population can have a sig- 
nificant impact on the inferred integrated SFR and the 



bright end of the IR LF. Although the fraction of ex- 
pected AGN in a cluster may be as low as a few percent, 
even a small number of AGN can dramatically increase 
the measured global SFR of a cluster, since AGN are 
preferentially IR bright compared to the normal star- 
forming population. We also warn against excluding red 
sequence galaxies in an attempt to remove AGN from 
a mid-IR sample of galaxies. In the Bullet Cluster, we 
have found that 20% of the non-AGN MIPS members lie 
on the red sequence. 

In addition, it is important when interpreting the in- 
tegrated specific SFRs of clusters that one has a good 
understanding of the infalling galaxy population. In the 
case of the Bullet Cluster, a few galaxies most likely be- 
longing to an infalling group heavily drive the integrated 
specific SFR, making the Bullet Cluster one of the more 
"active" clusters in Figure |H1 We have been able to iden- 
tify the important effect of the infalling population by 
comparing the IR LF of the Bullet Cluster to other sys- 
tems, and knowing the mass ratio of the infalling group to 
main cluster. In cases where the cluster dynamics are not 
well characterized, particularly in cluster mergers where 
infalling groups may be prevalent, it is important to con- 
sider that the integrated specific SFR is not necessarily 
reflective of effects from the merger itself, but may be 
dominated by an infalling group population. 

4.6. Transformation Mechanisms 
4.6.1. Timescale constraint from IR LF 

The IR LF of the Bullet Cluster reveals an excess of IR 
star-forming galaxies relative to other massive clusters. 
However, the observed Bullet Cluster IR LF is consistent 
with a combined cluster and group IR LF (appropriately 
scaled for mass) within a factor of ~1.5. This supports 
the hypothesis that galaxies from the infalling group pop- 
ulation (or "subcluster") are predominantly responsible 
for the enhanced star formation observed in the Bullet 
Cluster relative to other clusters. In this case, we can 
place a lower limit on the timescale required for a cluster 
mechanism to quench star formation in recently acquired 
group galaxies. Since subcluster star-forming galaxies 
are observed ~250 Myr after core passage, we can rule 
out fast-acting mechanisms such as ram pressure, as be- 
ing a dominant physical process in driving the evolution 
of star formation rate in a cluster merger. Instead we re- 
quire a slow-working process that can explain the excess 
of star-forming LIRGs observed in the Bullet Cluster 250 
Myr after a major merger. 

The need for a cluster process that works to quench 
star formation over long timescales of a few Gyr has been 
supported by much observational and theoretical work. 
There are two main physical mechanisms that work over 
such timescales in the cluster environment - galaxy ha- 
rassment and strangulation (also known as starvation). 
Of these two, we prefer strangulation to be the more 
plausible mechanism to eventually quench the remain- 
ing star formation in the Bullet Cluster and truncate the 
bright end of the IR LF. Strangulation is a slow process 
that begins to work in the cluster outskirts. The inter- 
action between the cluster ICM and the loosely bound 
hot halo gas of a galaxy can cause the halo to be gen- 
tly pushed out, eventually leading to transformati ons in 
both morphology and star formation ([Larson et al.lH980t 
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iBaloghet all [20001. Simulations by iBekki et all pOOl 
have shown that ram pressure and global tidal effects can 
be an effective means of removing halo gas within a few 
billion years. Once this halo gas is removed, the galaxy 
no longer has a source of replenishment for cold gas to 
continue star formation. 

There is much observational evidence in support of 
strangulation as an important cluster mechanism to ex- 
plain the dearth of late type star-f orming galaxies in 
cluster cores. This includes work by iTreu et all (2003) 
who classified morphologies of over 2000 galaxies in 
C10024+16 out to a 10 Mpc diameter region, and dis- 
covered "mild gradients" in the morphological fractions 
at large radii. They argue that such a trend could only be 
explai ned by a slow working mec hanism such as strangu- 
lation. iMoran et al.1 (J2006, 2007) also support strangula- 
tion as an important cluster process, based on the spatial 
distribution of quenched spiral galaxies and using ultra- 
violet and spectroscopic signatures of star formation to 
show that star formation in passive spirals must have 
decayed over ^S 1 Gyr. 

In addition to the observational studies, there is also 
mounting theor etical evidence in fav or of strangulation, 
such as work bv lBalogh et a l. (2000) who are able to re- 
produce the systematic observed difference in colors of 
cluster versus field galaxies, if they assume a gradual 
decline of SFR over a few Gyr in the galaxies accreted 
into the cluster. It has al s o bee n shown by simulations 
from lKawata &: Mulchaeyl i|2008l ) that strangulation can 
occur even in low mass groups, where the ICM-galaxy 
halo i nteraction is weaker than in clusters. More re- 
cently, iMcGee et al.l (|2009f ) have demonstrated through 
semi-analytic galaxy formation models that strangula- 
tion is a favorable mechanism to quench star formation 
over long timescales in order to match the observed evo- 
lution of red galaxy fractions in clusters. Further theo- 
ret ical evidence in favor o f strangulation includes work 
bv lWeinmann et"aLl (|2009T ) who have demonstrated that 
models including slow removal of hot halo gas in a galaxy 
can match observed fractions of passive cluster satellite 
and central galaxies. 

Galaxy harassment is also a mechanism that works 
on long timescales of a few Gyr from the cumulative 
effect of high-speed near encounters between galaxies, 
which can disturb galaxy morphologies and gas content 
(jMoore et al.lll996tlHavnes et aLll2007H . It is a process 
most efficient in cluster cores, where galaxy densities 
are high and near encounters occur frequently. How- 
ever, all of our LIRGs lie beyond R ~ 1 Mpc. Assum- 
ing that the LIRGs take ~1 Gyr to reach the cluster 
core where harassment can begin to have an impact, it is 
likely that other ICM processes effective in the interme- 
diate/outskirts region (such as strangulation or galaxy- 
galaxy interactions) would have already begun to take ef- 
fect. In addition, it has been suggested that harassment 
is only mildly effective in the case of dwarf cluster galax- 
ies (Mtot ~ 10 10 M^), with o nly 10% stellar mass loss 
on average ([Smith et al.1 120101 ) . Given that harassment 
has such a negligible effect on dwarf cluster galaxies, it 
is unlikely that it can significantly quench star formation 
in intermediate mass LIRGs. 

4.6.2. Spatial Distribution of LIRGs 



In addition to using timescale constraints based on the 
comparison of the Bullet Cluster IR LF to other systems, 
we can also look at the spatial distribution of the IR 
bright galaxies in the Bullet Cluster to constrain mecha- 
nisms that may play a key role in triggering or quenching 
star formation in cluster mergers. Although the Bullet 
Cluster shows a higher number density of IR galaxies 
compared to other clusters at all Lm values in Figure [TOl 
it is particularly useful to examine the spatial distribu- 
tion of the LIRGs/ULIRG. Looking at the spatial distri- 
bution of these LIRGs/ULIRG is a good way to isolate 
the galaxies that most likely do not belong to a "normal" 
galaxy cluster IR LF. This is evident in Figure [T01 which 
shows that the IR LFs of evolved Coma and CL1358+62 
truncate at a much lower Lju than the Bullet Cluster IR 
LF. 

The spatial distribution of the LIRGs and ULIRG, il- 
lustrated in Figure [T3J shows that all six IR bright galax- 
ies lie beyond the central Mpc region, where a Mach 3 
shock front propagates through the main cluster galax- 
ies, creating a strong ram pressure environment. There 
is also no clear correlation between the distribution of 
the LIRGs and the cluster merger structure as shown 
by X-ray surface brightness contours. This implies that 
the LIRGs were not triggered by the merger event itself, 
and in particular not trigge red by ram p r essure . This is 
consistent with results from lChung et all (J2009D who an- 
alyzed the specific SFRs of post-shock versus pre-shock 
galaxies in the Bullet Cluster and found that ram pres- 
sure did not significantly trigger or quenc h recent star 
forma tion in the Bullet Cluster galaxies. [Haines et al.l 
(2009) similarly concluded that cluster mergers do not 
have a strong impact on star-forming IR galaxies by 
finding no correlation between the dynamical state of a 
cluster (whether it is relaxed or merging), and the frac- 
tion of IR lu minous galax i es in a sample of 30 clusters. 
In contrast, IMoran et al.l (|2005) found evidence in fa- 
vor of triggered star formation traced by [Oil] emission 
in faint early type galaxies residing near the virial ra- 
dius of CL0024+16, where cluster merger effects may 
have genera ted shocks and/or induced galaxy harass- 
ment. While IMoran et all (|2005| ) favor a fast mechanism 
such as merger-induced ram pressure to trigger star for- 
mation in galaxies near the cluster outskirts, they also 
require a gradual fading of star formation on timescales 
of ~1 Gyr based on the radial distribution of Balmer ab- 
sorption line galaxies in CL0024+16. This is similar to 
the long timescale required to quench star formation in 
the Bullet Cluster LIRGs, if they are indeed from the 
subcluster population. 

Galaxy-galaxy mergers or tidal interactions also do not 
explain the LIRG population in the Bullet Cluster, de- 
spite the preferential distribution of LIRGs in the out- 
skirts, where the chance of mergers is relatively high 
within a cluster environment. We find no obvious signs of 
recent interactions that would have triggered star forma- 
tion from th e morp hologies of the LIRGs and ULIRG. 
IGeach et al.l (|2009f ) found that in the cluster merger 
CL0024+16, several LIRGs belong to smaller groups that 
may be falling into the cluster, with some evidence that 
star formation is triggered by interactions within these 
small groups. However, they also find no statistical dif- 
ference in the local environment of LIRGs versus quies- 
cent spiral galaxies. 
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The current data thus provide no evidence that the ob- 
served LIRGS are triggered by either the cluster merger 
or galaxy-galaxy interactions. It also seems implausible 
that these LIRGS are field galaxies falling into the clus- 
ter for the first time, based on a comparison of the IR 
LFs of the Bullet Cluster with those of evolved Coma 
and CL1358+62. The IR LFs of both evolved Coma and 
CL1358+62 truncate at lower values of Lju than in the 
Bullet Cluster IR LF, which would be unlikely if the Bul- 
let Cluster LIRGs were inf ailing field galaxies. 

Figure [T3] shows that the LIRGs are all distributed 
at R > 1 Mpc, consistent with recent studies which 
find that mid-IR and sub- mm dusty star-forming galaxies 
avoid the cluster core and are preferentially distributed in 
the in termedia te and outsk irts regions fe.g. lBraglia et all 
l2010t : Kovama et al.l 120101 ) . However, the spatial distri- 
bution of the LIRGs shown in Figure H"2l also corresponds 
to the outskirts region of the subcluster, where it is more 
likely to find dusty star-forming galaxies. All the LIRGs 
are within 1 to 2 times the -R200 radius of the subcluster, 
from the center of the "bullet" . 

Overall, the spatial distribution of the Bullet Cluster 
LIRGs shows no evidence of being strongly influenced by 
the cluster merger. Combined with results from compar- 
ing the IR LF of the Bullet Cluster to the IR LFs of 
other systems, the LIRGs are unlikely to be field galax- 
ies. However, the spatial distribution of the LIRGs is 
consistent with the galaxies belonging to the outskirts 
region of the subcluster. 

5. CONCLUSIONS 

We use 24/im MIPS data to calculate the mass nor- 
malized global star formation rate of 37 spectroscopi- 
cally confirmed Bullet Cluster members within R < 1.7 
Mpc, excluding AGN sources using IRAC, X-ray, and op- 
tical spectroscopic data. The integrated obscured SFR 



of the Bullet Cluster is 267 M 



©yr 



with a 0.2 dex un- 



certainty, using the L24-SFR relation from iRieke et al.1 
(2009) . Normalizing by cluster mass we obtain a specific 
star formation rate o f 28 M Q yr _1 per 1 14 M(ry Us i ng the 
iDale fc Heloul (j2002f ) models with the iKennicuttl (|1998f ) 
SFR calibration, and applying a lower limit on 24/im flux 
comparable to previous works, we obtain a mass normal- 
ized integrated SFR of 47 M @ yr -1 per 10 14 M Q . We find 
that the specific SFR of the Bullet Cluster is one of the 
highest among other merging and non-merging clusters 
from the literature, with the exception of CL0024+16. 
In the Bullet Cluster, a small population of LIRGs and 
ULIRG contribute 30% and 40% of the total obscured 
star formation rate. These galaxies comprise the bright 
end of the infrared luminosity function, which are not 
observed in the IR LFs of evolved Coma or CL1358+62. 
A Schechter fit to the IR LF of the Bullet Cluster yields 
L* = 44.70±0.16, which is -0.25 dex brighter than L* of 
evolved Coma and 0.35 dex brighter than CL1358+62, a 
massive cluster at z = 0.3. We observe an excess of 
IR bright sources in the Bullet Cluster along the entire 
range of infrared luminosities in the IR LF, compared to 
the IR LFs of evolved Coma and CL1358+62. One way 
to explain this is to attribute the excess of IR sources 
to the galaxy population associated with the group that 
has recently been accreted into the main cluster —250 
Myr ago (core passage) . The combined IR LF of evolved 
Coma and mass-scaled SG1120 is consistent within a fac- 



TABLE 1 
Data 






Name 


MIPS 


MIPS AGN 


ALL 


Confirmed members 

Interlopers 
Candidate members 


44 
83 
14 


6 

5 
2 


406 
495 



Note. — The number of galaxies that are spectro- 
scopicaly confirmed cluster members, interlopers and 
cluster candidates (based on optical and IRAC colors). 
The columns show number of galaxies from the MIPS 
catalog with WFI and IRAC counterparts, number of 
AGN within the MIPS catalog determined by IRAC col- 
ors, X-ray emission, and optical spectroscopy, and the 
total number of galaxies from the entire IMACS spec- 
troscopic catalog. 
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Fig. 1. — Spectroscopic redshift distribution with shading to indi- 
cate cluster member region. Inset shows zoomed-in view of cluster 
members as defined by the caustic analysis of the infall region. 

tor of -1.5 of the observed Bullet Cluster IR LF. This 
supports the idea that the galaxies originating from the 
group population have not yet been processed to match 
typical cluster levels of star formation rate. In this case, a 
transformation mechanism such as strangulation, which 
acts on timescales longer than the merger timescale of the 
main cluster and subcluster ( ^ 250 Myr) is necessary to 
eventually truncate the bright end and suppress the faint 
end of the IR LF to match the observed distribution in 
Coma and CL1358+62. 
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JPL/Caltcch. The authors acknowledge support for this 
work from NASA/Spitzer grants 1319141 and 1376614. 

Facilities: Spitzer (IRAC), Spitzer (MIPS), CXO 
(ACIS-I), MagellamBaade (IMACS) 



Infrared Luminosity Function and Star Formation Rate of the Bullet Cluster 



11 





TABLE 2 








LIRGs 


AND AGN 






Object Type 


RA [hh:mm:ss] 


Dec [° ' "] 


l ir [L ] 


SFR [M yr- 1 ] 


ULIRG/LINER 


6:58:30.87 


-56:03:36.31 


1.3 x 10 12 


104 


LIRG 


6:59:13.45 


-55:56:31.19 


2.5 x 10 11 


28 


LIRG 


6:58:55.33 


-55:55:43.04 


1.5 x 10 11 


16 


LIRG 


6:58:29.30 


-55:53:15.92 


1.4 x 10 11 


11 


LIRG 


6:57:47.86 


-55:55:10.83 


1.1 x 10 11 


11 


LIRG 


6:58:37.06 


-56:00:45.66 


1.0 x 10 11 


10 


LIRG, X-ray AGN 


6:58:03.52 


-56:01:14.23 


1.0 x 10 11 




X-ray & IRAC AGN 


6:58:35.22 


-56:01:4.78 


8.7 x 10 10 




X-ray AGN 


6:58:17.48 


-56:02:49.35 


8.1 x 10 10 




Optical AGN (Seyfert) 


6:58:42.50 


-56:00:28.22 


7.6 x 10 10 




IRAC AGN 


6:58:08.48 


-55:53:36.45 


1.7 x 10 10 




X-ray AGN 


6:58:26.67 


-56:00:0.18 


9.2 x 10 9 





Note. — Spcctroscopically confirmed members with MIPS emission that are cither LIRGs 
and/or AGN. Column 1 notes whether a galaxy is a LIRG/ULIRG, and if it is detected as 
cither an X-ray point source or an IRAC AGN wedge source. Columns 2-5 show the right 
ascension, declination, total infrared luminosity, and star formation rate for each galaxy. SFR 
of the ULIRG assumes that 35% of the total IR luminosity is powered by AGN, and the rest 
by star formation. For the LIRGs and AGN, the lcr error of the measured 24/xm flux is 
between 1-10%, and thus the final error of the total infrared luminosity is dominated by the 
uncertainty in the L24 to Lir calibration, which is ^S 0.15 dcx. 
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Fig. 2. — The R-[4.5] and B-R colors of spectroscopically confirmed cluster members (dots), cluster members with MIPS emission (solid 
circle), cluster candidates with MIPS emission (open circle), and known interlopers with MIPS emission (cross). The spectroscopically 
confirmed cluster galaxies, particularly those with MIPS emission, form a tight sequence in this color-color space. To minimize contamination 
from background MIPS sources, we include galaxies only within the diagonal boundaries indicated. Galaxies with B — R > 3 are also 
excluded since these are optically redder than the red sequence. The right panel shows the R-[4.5] and B-R colors of confirmed cluster 
members only. Five confirmed cluster LIRGs are shown as blue stars, and one ULIRG as a green star with a square border. 
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Fig. 3. — IRAC color-color diagram, showing the AGN wedge (Lacy et al. 2004; Stern ct al. 2005). Symbols are the same as for Figure |21 
MIPS sources identified as IRAC or X-ray AGN are shown with a blue square and plus sign overplotted on a solid or open red circle. The 
Seyfert galaxy identified in Figure [4] is shown with a blue diamond behind the ULIRG (green star symbol). The dotted line shows the M82 
track starting from z = then progressing to z = 0.3,0.5, 1, and 2, with each redshift marked by a cross. It is clear that many of the high 
redshift background sources are removed from our color selection based on Figure [2] 
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TABLE 3 

SCHECHTER PARAMETERS 



Cluster/group 



logL* 4>* [N/Mpc 2 /log L IR ] 



Coma (evolved z = 


0.3) 


-1.4 


44 44+0-27 
**•**— 0.24 






CL1358+62 




-1.4 


44.331°,;^ 


xm+\% 




SG1120 




-1.4 


44.991°,;" 


i.25t°; 5 3 




Bullet Cluster 




-1.4 


44.681°,- 


-, 77 +0.5 


1 

CD 


Bullet Cluster (with AGN) 


-1.4 


44.8iig;i| 


1 3fi+ - 4 





NOTE. — Schcchtcr parameters for the clusters and S G1120 shown in Fig- 
urc llOl Schcchter parameters for Coma arc adopted from Bai cl al. (2009) and 
evolved to z — 0.3 usi ng L* oc (1 + z) e volution, while values for CL1358+62 
and SG1120 are from ITran et all J2009D . We fit a Schcchter function to the 
observed Bullet Cluster IR LF, without and with the AGN, fixing the faint 
end slope a to the match that of Coma. 
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one ULIRG are highlighted as star symbols, with a square border 
around the ULIRG. 




Fig. 7.— The ULIRG imaged with HST/ACS in the F606W filter. 
It is a barred spiral with no signs of recent major interactions. 



Fig. 8. — The optical color-magnitude diagram with spectroscopi- 
cally confirmed cluster members (dot), spectroscopically confirmed 
members with MIPS emission (solid circle) and cluster candidates 
with MIPS emission (open circle). MIPS sources identified as AGN 
with IRAC, X-ray, or optical data are indicated with a blue square, 
triangle, and diamond with a plus sign, respectively. Five spec- 
troscopically confirmed LIRGs are shown as blue stars, and one 
ULIRG/LINER as a green star with a square border. Grey solid 
line shows the col or-magnitude relation wit h slope and normaliza- 
tion adopted from Lopez-Cruz ct al. (2004). 
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Fig. 9. — The mass normalized integrated SFR as a function 
of redshift, with th e Bullet Clust e r sho wn as a blue solid circle. 
SFRs calculated by IHaines et"aH 1 20091) for A1758 is illu s trated 
with a triangle. All others are calculated in Geach ct al.. (120061) 
and references therein. 



Infrared Luminosity Function and Star Formation Rate of the Bullet Cluster 



15 



1.5 



1.0 



ST 0.5 



0.0 



-0.5 



-1.0 



»'■.-' ' 



"Y* 



AGN Excluded 



** 



\4. 



• Bullet Cluster 
ACL1358 (Tron + 09) 
5KSG1120 (Tran + 09j 
iComo ot z = 0.3 (Boi- 



09) 



1.5 



1.0 



<T 0.5 



0.0 



-0.5 



-1.0 



TT7r7 



I ' ' :' I 

:AGN Not Excluded 






H 



• Bullet Cluster 
iCL1358 (Tran + 09) 
*SG1120 (Tran + 09J 
■ Coma ot z = 0.3 (Bai+09) 



\ * 



42 



43 44 45 

Log L,„ [ergs/s] 



46 



42 



43 44 45 

Log L,„ [ergs/s] 



46 



Fig. 10. — The infrared luminosity function of the Bullet Cluster (solid circle), Coma cluster (solid square), CL1358+62 (solid triangle), 
and SGI 120 (asteris k), with the IR LF of the Coma cluster evolved to z ~ 0.3. All infrared luminosities shown here are based on the 
I Dale & Heloul H2002T ) galaxy templates. The Schechter fit to the Bullet Cluster IR LF is shown as the dashed curve, and the MIPS 80% 
completeness limit is indicated as a dotted vertical line. The Left panel shows the Bullet Cluster IR LF excluding all known AGN candidates 
and the right panel includes six AGN confirmed members and two AGN candidate members (scaled by 0.35 to account for probability of 
being actual cluster members). The inclusion of just a small number of AGN can significantly elevate the IR LF, preferentially in the bright 
end. However, we note that even without the AGN, the Bullet Cluster IR LF exhibits an excess of IR luminous sources relative to Coma 
and CL1358+62. The brightest Bullet Cluster galaxy is a ULIRG and is not included in the Schechter fit. 
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Fig. 11. — The infrared luminosity function of the Bullet Cluster 
(solid circle) excluding known AGN. Overplotted with solid (red) 
square is the IR LF of Coma evolved to z = 0.3 added with the 
IR LF of SG1120, after the latter has been scaled down by a fac- 
tor of ~4 to match the approximate mass of the infalling group 
population in the Bullet Cluster. 




Fig. 12. — R-band WFI image shown with X-ray surface bright- 
ness contours. The large box shows the FOV of our 24/^m MIPS 
data, and small boxes indicate the location of the five LIRGs and 
one ULIRG. 
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